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Mitigating  Error  Propagation  Effects  in  a  Decision 

Feedback  Equalizer 

Michael  Reuter,  Member,  IEEE ,  Jeffery  C.  Allen,  Member,  IEEE ,  James  R.  Zeidler,  Fellow,  IEEE ,  and 

Richard  C.  North,  Member,  IEEE 


Abstract — We  present  an  approximate  analysis  approach  to  the 
computation  of  probability  of  error  and  mean  burst  error  length  for 
a  decision  feedback  equalizer  (DFE)  that  takes  into  account  feed¬ 
back  of  decision  errors.  The  method  uses  a  reduced-state  Markov 
model  of  the  feedback  process  and  is  applicable  to  linear  modula¬ 
tion  formats.  We  use  this  technique  to  analyze  a  DFE  design  that 
mitigates  the  effects  of  feedback  error  by  incorporating  a  soft  de¬ 
cision  device  into  the  feedback  path  and  a  norm  constraint  on  the 
feedback  filter  weights.  We  apply  the  DFE  design  and  analysis  ap¬ 
proach  to  a  dispersive  multipath  propagation  environment. 

Index  Terms — Decision  feedback  equalizers,  error  analysis, 
error  propagation,  Markov  processes. 

1.  INTRODUCTION 

HE  DECISION  feedback  equalizer  (DFE)  is  an  impor¬ 
tant  component  in  many  digital  communication  receivers 
and  is  used  to  suppress  intersymbol  interference  (ISI)  eaused 
by  dispersive  propagation  channels  [1],  [2],  as  well  as  rejeet 
in-band  interference  [3],  [4].  The  DFE  incorporates  a  feedfor¬ 
ward  filter  that  operates  on  the  received  signal  to  suppress  pre¬ 
cursor  ISI,  with  a  feedback  filter  that  operates  on  previously 
detected  ehannel  symbols  to  suppress  posteursor  ISI.  The  DFE 
generally  outperforms  the  traditional  linear  equalizer,  particu¬ 
larly  if  the  ehannel  has  deep  spectral  nulls  in  its  response  [1] 

However,  degradation  in  DFE  performance  occurs  when  m- 
eorreetly  detected  symbols  are  fed  through  the  feedback  filter. 
Then  instead  of  mitigating  ISI  from  the  eursor  sample,  the  DFE 
enhances  ISI.  Error  propagation  may  result  that  eauscs  bursts 
of  decision  errors  and  a  corresponding  increase  in  the  average 
probability  of  bit  and  symbol  error.  Moreover,  the  bursty  nature 
of  DFE  errors  has  implications  for  error  correction  coding  and 
interleaver  depth  that  may  be  incorporated  into  the  receiver  de- 
sign  [5],  [6], 

Error  propagation  generally  does  not  severely  affeet  DFE  per¬ 
formance  if  the  ehannel  delay  spread  is  on  the  order  of  a  symbol 
duration  or  less  [  1  ],  [7].  However,  as  bandwidths  increase  to  ac¬ 
commodate  high  data  rates,  the  associated  increase  in  relative 
delay  spread  may  result  in  error  propagation  for  a  given  ehannel. 
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Also,  high-order  modulation  formats  may  be  more  suseeptible 
than  low-order  formats.  Therefore,  a  need  exists  to  quantify 
DFE  performance  in  the  presence  of  decision  error  feedback  and 
to  mitigate  the  effects  of  error  propagation 

Performance  analysis  of  a  DFE  that  includes  decision  error 
feedback  is  a  difficult  problem  and  has  received  considerable 
theoretical  attention.  Using  an  independence  assumption  on  the 
noise  and  assuming  no  precursor  ISI,  Monsen  used  a  Markov 
model  of  the  error  feedback  proeess  to  derive  probability  of  error 
for  the  DFE  [1].  The  mean  burst  error  length  also  ean  be  de¬ 
rived  using  this  theory  [8].  Although  this  approach  is  accurate, 
it  becomes  computationally  burdensome  beeause  the  number  of 
Markov  states  exponentially  increases  as  the  number  of  feed¬ 
back  taps  increases.  Duttweiler,  Mazo,  and  Messersehmitt  sim¬ 
plified  the  analysis  by  aggregating  the  Markov  states  into  a  re- 
dueed-state  model  and  then  bounding  the  state  transition  proba¬ 
bilities  [9].  Bounds  on  the  probability  of  error  [9]  and  the  burst 
error  length  [10]  are  obtained  Tighter  bounds  are  possible  by 
refining  the  definition  of  the  aggregate  states  [1 1],  [12]. 

Various  techniques  for  mitigating  error  propagation  also  have 
been  proposed.  Large  taps  in  the  feedback  filter  of  an  equalizer 
designed  using  the  minimum  mean-square-error  (MSE)  crite¬ 
rion  cause  significant  self-generated  ISI  if  decision  errors  are 
made.  A  potential  solution  is  to  jointly  optimize  the  feedfor¬ 
ward  and  feedback  filters  with  an  additional  constraint  on  the 
norm  of  the  feedback  filter  taps  This  regularization  approach 
was  applied  to  the  magnetic  recording  ehannel  in  [13],  [14]. 
DFE  structures  also  have  been  proposed  recently  that  eontain 
a  soft  decision  device  in  the  feedback  path  to  compensate  for 
unreliable  decisions  [15] — [  1 8]. 

Combining  these  techniques,  we  present  a  minimum  MSE 
DFE  design  that  incorporates  both  a  discrete  soft  decision  device 
and  a  norm  constraint  on  the  feedback  filter  to  mitigate  error 
propagation.  Reliable  deviee  inputs  that  are  close  with  respect 
to  the  L 2  norm  to  a  component  of  the  symbol  constellation  are 
fed  back  as  hard  decisions.  Unreliable  deviee  inputs  that  are 
not  close  are  fed  back  as  intermediate  decisions.  Also,  feedback 
gain  is  controlled  by  the  norm  constraint. 

To  capture  the  bursty  nature  of  the  DFE  errors,  we  use  both 
the  probability  of  symbol  error  P e  and  the  mean  burst  error 
length  fiBL  as  the  primary  DFE  performance  measures.  We  de¬ 
termine  these  quantities  analytically  using  the  fundamental  ap¬ 
proach  of  [  1  ],  but  we  derive  an  approximate  Markov  model  with 
significantly  fewer  states  by  making  an  independence  assump¬ 
tion  involving  the  feedback  error  sequence.  We  demonstrate 
that  design  trade-offs  between  Pe  and  /xbl  result  in  a  classic 
“L-eurve.”  From  this  geometric  interpretation,  an  optimum  de¬ 
sign  becomes  evident. 
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Fig.  1  Receiver  front  end. 

II.  System  definitions 

We  are  concerned  with  coherent  demodulation  of  data  trans¬ 
mitted  through  a  discrete  multipath  environment.  The  channel 
symbol  modulation  is  M- ary  phase  shift  keying  (A^-PSK).  We 
operate  on  the  complex  baseband  signal  using  the  receiver  front 
end  shown  in  Fig.  1 .  The  down-converter  is  locked  to  the  car¬ 
rier  frequency  fc  but  not  the  phase.  Therefore,  any  phase  com¬ 
pensation  is  done  by  the  equalizer  [19].  The  noise  at  the  input 
to  the  ideal  lowpass  filter  is  modeled  as  a  circularly  symmetric 
complex  white  Gaussian  process  with  two-sided  power  spectral 
density  N q. 

The  signal  at  the  output  of  the  ideal  lowpass  filter  is  modeled 

as 

L  oo 

x(t)  =  Pi(t)  ^2  s{m)Pr{t  -  thT  -  Aj  -  r) 

1=1  m=— oo 

xe-<2  W/.(A,+r)  +  n(t)  (1) 

where 

L  number  of  discrete  multipath  components; 

A i  relative  path  delays; 

s(m)  channel  symbols; 

T  symbol  duration. 

The  multiplicative  fading  terms  fix (t)  incorporate  all  channel 
effects  and  are  modeled  as  mutually  independent  random 
processes  with  zero  mean  and  power  The  corresponding 
path  delays  A i  are  deterministic.  The  symbol  pulse  pr{t)  is 
a  square-root  raised  cosine  waveform  with  unit  energy.  The 
variable  r  represents  symbol  timing  error  and  is  modeled  as 
a  continuous  random  variable  uniformly  distributed  between 
—T/2  and  T/2.  The  A/-PSK  symbols  s(m)  are  mutually 
independent  random  variables  with  unit  power.  The  noise 
n(£)  is  a  circularly  symmetric  Gaussian  process  with  power 
<j\  —  2 Nq/T.  Using  these  assumptions  and  by  averaging  over 
all  random  quantities,  we  define  the  received  signal-to-noise 
ratio  (SNR)  to  be  [20] 

^NR=2iVo^^''  (2) 

The  fading  on  all  paths  is  slow  with  respect  to  T  so  that  we  can 
assume  a  quasistatic  channel  with  a  given  delay  structure.  All 
performance  results  are  based  upon  the  minimum  MSE  Wiener 
filter  that  is  optimized  with  respect  to  a  static  realization  of  the 
channel  defined  by  the  set  of  parameters  {/?i,  fa,  •  •  • ,  fth  r}- 
Also,  the  bandwidth  of  the  lowpass  filter  in  Fig.  I  accommo¬ 
dates  sampling  x(t)  at  the  rate  2/T  to  get  the  discrete  samples 
x(kT/ 2).  This  allows  us  to  incorporate  a  fractionally  spaced 


Fig.  2.  DFE  structure  with  soft  decision  feedback  and  hard  decision  output. 

feedforward  filter  into  the  DFE  design  to  provide  robustness 
against  timing  error  r  as  well  as  matched  filtering  to  the  trans¬ 
mitted  pulse  pr{t)  and  the  propagation  channel  [2]. 

111.  DFE  Design 

A  schematic  of  the  DFE  design  is  given  in  Fig.  2.  The  sym¬ 
metric  fractionally  spaced  feedforward  filter  wp  has  2 N  +  1 
taps  with  N  taps  leading  and  N  taps  lagging  the  cursor  tap  The 
symbol-spaced  feedback  filter  wb  operating  on  the  output  of  the 
soft  decision  device  has  N&  taps.  We  define  the  feedforward  and 
feedback  data  vectors  as 

x(fc)  =  [x(kT  +  NT/2)  x(kT  +  (N  -  l)T/2) 

■  ■  ■  x(lcT)  •  •  •  x{kT  -  NT/2)Y 
d{k)  =  [d{k  -  1)  d(k  —  2)  •••  d{k-Na)]i  (3) 

where  f  denotes  the  matrix  transpose.  Then  the  input  to  the 
hard  decision  device  used  to  estimate  the  cursor  symbol  s(k) 
is  written  as 

y(k)  =  w  px(k)  +  wg  d  (fc)  (4) 

where  H  denotes  the  Hermitian  transpose 

There  are  many  ways  to  define  the  soft  decision  device  in 
Fig.  2.  We  propose  the  discrete  implementation  shown  in  Fig.  3 
for4-PSK.  Decision  regions  for  higher-order  Af-PSK  scale  ac¬ 
cordingly.  The  shaded  areas  represent  the  intermediate  decision 
regions  and  arc  defined  by  the  angle  0.  The  intermediate  deci¬ 
sions  are  averages  of  the  adjacent  complex  Af-PSK  symbols. 
The  intent  is  to  minimize  the  distance  between  the  intermediate 
decision  and  the  adjacent  symbols.  Applying  this  reasoning  to 
2-PSK  modulation  results  in  the  erasure  scheme  presented  in 
[16] 

We  derive  the  optimum  constrained  Wiener  weights  for  the 
DFE  by  assuming  error-free  past  decisions  obtained  from  a 
training  sequence  so  that  d(k)  =  s(k).  The  Wiener  realization 
of  the  DFE  is  the  solution  of  the  regularization  problem  [21  ] 

min{£[|s(fc)  -y(fc)|2]} 

subject  to  ||wD||2  =  7||wDu||2  (5) 

where  E[  ]  is  the  expectation  operator  taken  over  the  informa¬ 
tion  and  noise  processes,  and  ||wDu  ||  is  the  norm  of  the  uncon¬ 
strained  feedback  weights.  The  parameter  7  controls  the  norm 
constraint  where  0  <  7  <  1. 
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Fig.  3.  Soft  decision  regions  for  4-PSK  wilh  inicrmcdiaic  decisions  in  shaded 
regions. 

Using  the  method  of  Lagrange  multipliers,  the  optimum 
Wiener  weights  are  readily  shown  to  be  the  solution  of  the  set 
of  linear  equations  [2] 


R 

Q" 

wp  ' 

p 

.Q 

(1  +  A)I;vd  . 

.  wB. 

0 

> 

a 

where  R  =  £?[x(fc)x(fc)"]f  Q  =  E[d(k)x(k)H], 

p  =  E[$*  (k)x(k)]>  and  *  denotes  complex  conjugation. 
I/vB  is  the  identity  matrix  of  size  Nb,  0,vb  is  the  all-zero 
column  vector  of  length  iVo,  and  A  is  the  Lagrange  multiplier 
A  >  0.  Equation  (6)  shows  that  the  norm  constraint  analytically 
introduces  synthetic  white  noise  with  power  A  into  the  feedback 
path. 

The  optimum  feedforward  and  feedback  filter  weights  are 
written  explicitly  as 

wF=(R-aQ"Q)-lp 
wD=-aQwr  (7) 

where  a  =  1/(1  4-  A)  and  0  <  a  <  1.  If  a  =  1,  the  optimum 
weights  represent  the  standard  unconstrained  DFE.  As  a  — » 
0,  the  feedback  path  diminishes  and  the  feedforward  filter  wF 
approaches  the  optimum  linear  equalizer  [13]. 

For  convenience,  we  will  define  the  norm  constraint  in  terms 
of  the  parameter  7  in  (5).  However,  there  is  no  closed-form  re¬ 
lationship  between  a  and  7  except  at  the  extreme  values  where 
7  =  0  and  7  =  1  correspond  to  a  =  0  and  a  =  1,  respectively. 
Therefore,  when  invoking  a  constraint  with  arbitrary  7,  we  use 
(7)  in  a  recursive  search  strategy  to  find  the  corresponding  a. 
Although  we  are  not  concerned  with  implementation  issues,  we 
note  that  if  the  DFE  weights  are  adapted  with  the  least  mean 
square  algorithm  (LMS),  the  constrained  equalizer  can  be  im¬ 
plemented  with  only  a  minor  increase  in  complexity  using  leaky 
LMS  [13],  [22,  pp.  746-747] 


IV.  DFE  Performance  Analysis 

We  use  the  average  probability  of  symbol  error  PE  and  mean 
burst  error  length  as  the  two  primary  measures  of  perfor¬ 
mance  for  the  DFE.  The  approach  we  take  to  compute  these  two 
values  is  to  separate  the  decision  device  input  y(k)  into  com¬ 
ponents  representing  the  self-generated  ISI  resulting  from  de¬ 
cision  error  feedback,  the  remaining  residual  ISI  not  spanned 
by  the  feedback  filter  wg,  and  noise.  We  further  separate  the 
error  feedback  sequence  of  length  TVq  into  the  actual  transmitted 
sequence  and  a  multiplicative  process  that  represents  potential 
A/-PSK  phase  shift  errors  at  the  output  of  the  soft  decision  de¬ 
vice.  Wc  then  make  a  simplifying  assumption  that  this  multi¬ 
plicative  error  process,  the  transmitted  symbols,  and  the  noise 
process  are  statistically  independent  of  each  other.  An  approxi¬ 
mate  Markov  model  of  the  DFE  follows  in  which  each  Markov 
state  is  defined  by  a  permutation  of  the  multiplicative  error  se¬ 
quence.  Using  this  independence  assumption,  we  calculate  com¬ 
ponents  of  the  corresponding  Markov  probability  transition  ma¬ 
trix  by  averaging  over  the  noise  and  all  possible  permutations  of 
the  transmitted  A/-PSK  symbols  contained  within  the  sample 
y(k)  using  a  computational  technique  presented  in  [23].  Then 
straightforward  computations  produce  the  approximate  values 
PE  and  /ibl- 

The  analysis  begins  by  decomposing  the  feedforward  data 
vector  using  (1)  and  (3)  as 

L 

x(fc)  =  ^2  ATis(fc)  4-  n  (k)  (8) 

;=i 

where  s (k)  is  the  infinite-length  symbol  vector 

s(k)  =  [•  •  •  s(k  +1)  s(k)  s(k  —  1)  •  •  •]*  (9) 

and  n (k)  is  a  (2N  -f  1)  x  1  white  complex  Gaussian  vector 
process  with  power  <j\.  Each  column  of  the  (2JV  4-  1)  x  00 
channel  matrix  T/  is  a  replica  of  the  fractionally  sampled  wave¬ 
form  pr(t  —  A  1  —  r)  exp(—i27tfc(Ai  4-  r)),  down-shifted  by 
two  samples  with  respect  to  the  column  on  the  left  and  truncated 
by  the  observation  window  of  length  2N  +  1.  By  construction, 
the  equalizer  is  eoarse-synehronized  to  the  path  with  A;  =  0  so 
that  the  center  row  of  each  T/  is  given  by 

Tl(N  +  1,  :) 

=  [•••  Pt(T  +  Aj  +  t)  Pt(Ai  +  t) 

Pt(—T  +  Aj  +  t)  . .  ,]e-2*/c<^+-0,  (io) 

Because  the  computational  technique  in  [23]  requires  a  finite 
number  of  ISI  symbols,  we  truncate  each  T /  to  get  the  (2N  4- 
1)  x  (2cZ  4  1)  matrix  T/  such  that  wc  maintain  d  columns  on 
either  side  of  the  column  containing  the  sample  Pt(^i  +  T) 
corresponding  to  the  cursor  symbol  s(k)  in  ( 1 0).  Because  of  the 
fast  roll-off  of  the  square-root  raised  cosine  pulse,  the  choice 
of  the  truncation  parameter  d  is  not  critical,  although  it  must 
be  chosen  with  respect  to  the  largest  multipath  delay  such  that 
d  >  max{|Aj|}/2’  and  d  >  iVB  (d  =  30  in  all  numerical 
results  presented  here). 
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Then  by  defining  the  composite  channel  matrix  T  as 


L 

T  =  X>T,  (11) 

/=  1 

the  feedforward  data  vector  is  well  approximated  by 

x(Ar)  =  Ts(Ar)  4-  n(fc),  (12) 

where  s (A;)  is  the  truncated  symbol  veetor  given  by 

s  (k)  =  [s(k  4-  d)  s(k  4-  d  —  1) 

•••  s(k)  •••  s{k-d)]l  (13) 


We  next  partition  T  and  s (k)  such  that 

'si(k) 


Ts  (Jfc) 

=  [Ei 

f 

=  Ei  si  (A;) 

where 

Si(k) 

=  [s{k 

+  d) 

s(k 

s2(k) 

=  {s(k 

-1) 

s(k 

s  3(k) 

=  [s(k 

-N0 

-1) 

E2  E3] 


s(k) 
S2  (k) 

Ls3(/c)  J 


s{k 


)  •  •  •  s(k  +  1)]* 

••  s(k  —  A'b)]1, 

Nb-2)  ■■■  s(k-d)]j. 

(15) 


Then  the  eross-eorrelation  components  comprising  the  optimum 
Wiener  weights  in  (7)  are  given  by  R  =  TTH  4- 
Q  =  E^,  and  P  =  f.  Also,  using  (4),  (7),  (12),  and  (14),  we 
ean  write  the  decision  deviee  input  as 

y(k)  =w f(x(&)  -  aQud(k)) 

=  Wp  (Ts(A;)  -  aE2d(A;))  4-  n(k) 

=  Wp  (E1S1  (A;)  +  fs(k)  4-  E3S3(A;)) 

+  WpE2(s2(A)  -  ad  (A:))  4-  Wp  n(fc).  (16) 

The  feedback  soft  decisions  ean  be  represented  as 

d(k)  =  r(k)Os2(k)  (17) 

where  0  denotes  element-wise  product  and  r(Ar)  is  a  random 
A^b  x  1  vector  in  which  each  term  represents  the  rotation  and 
sealing  of  the  output  of  the  soft  decision  deviee  with  respect  to 
the  transmitted  M -PSK  symbol  (Fig.  3).  If  there  are  no  decision 
errors,  r(  A:)  is  a  vector  of  ones  (r(  A:)  =  1a-d).  From  (16)  we  see 
that  if  the  norm  constraint  is  used  (a  <  1),  and  r(Ar)  =  l/vD, 
the  feedback  filter  does  not  cancel  all  of  the  postcursor  symbols 
within  its  span. 

The  feedback  error  term  in  (16)  then  can  be  written  as 

s2(A;)  -  ad(A;)  =  (1ab  -  ar(A;))  0  s2{k) 

=  e(A;)  0  s2(Aj)  (18) 

where  e(Ar)  is  the  random  multiplicative  error  vector.  The  pos¬ 
sible  values  for  each  element  of  e(A;)  are  displayed  in  Fig.  4  for 


Q 


I  e?  (i.o) 

x 
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Fig.  4  Multiplicative  feedback  error  values  for  4-PSK  with  soft  decisions 
(0  >  0)  and  a  norm  constraint  («  <  1). 

4-PSK  using  the  soft  decisions  (0  >  0)  given  in  Fig.  3.  If  0  =  0, 
the  intermediate  (even)  terms  in  Fig.  4  do  not  occur.  The  deci¬ 
sion  device  input  then  can  be  decomposed  in  terms  of  the  eursor 
symbol  ISI  and  noise  as 

y(k)  =  ao-s(A;)  4-  &iSi(k)  +  a2(e(A;)  0  s2(A;)) 

+a3s3(A;)  4*  Wp  n(Ar)  (19) 

where  no  =  Wp  f ,  and  the  remaining  ISI  coefficient  row  vectors 
are  given  by  a j  =  E j . 

The  approach  taken  by  Monsen  [1]  assumes  the  residual  ISI 
involving  the  terms  si(A:)  and  s3(A:)  in  (19)  is  negligible  Then, 
making  the  reasonable  approximation  that  the  noise  term  is 
white  Gaussian  and  therefore  independent  of  the  feedback  error 
proeess,  a  Markov  model  of  the  DFE  results  in  whieh  each 
state  is  defined  by  a  realization  of  the  feedback  error  sequence 
e(A;)  0  s2(A;).  By  conditioning  y(k)  on  the  cursor  symbol  s(k) 
and  each  realization  of  e(A;)  0  s2(At),  a  Gaussian  variate  results 
and  is  used  to  eompute  components  of  the  transition  probability 
matrix  P.  Although  wc  could  take  this  approach  by  expanding 
the  definition  of  each  state  to  include  realizations  of  the  residual 
ISI,  the  number  of  states  would  be  M2d~NB  {2M2)Nb  for 
A/-PSK  with  soft  decisions,  and  therefore  computationally 
unwieldy. 

The  approach  proposed  in  [24]  accomplishes  a  significant  re¬ 
duction  in  states  by  assuming  the  residual  ISI  not  spanned  by 
the  feedback  filter  is  statistically  independent  of  e(k)  0  s2(k) 
and  the  noise  The  states  are  defined  as  in  [1]  but  the  transi¬ 
tion  probabilities  are  computed  by  averaging  over  the  residual 
ISI  using  an  efficient  numerical  technique.  Adapting  this  ap¬ 
proximate  method  to  A/-PSK  with  soft  decisions  would  result 
in  (2A/2)Ad  states. 

We  achieve  a  further  reduction  in  states  (RIS)  by  making  the 
following  assumptions: 

RIS-1.  The  multiplicative  error  vector  e(A;)  is  statistically 
independent  of  the  noise  n(A;). 

RIS-2.  e(Ar)  is  independent  of  s(k)  and  the  residual  ISI 
represented  by  sa(A:)  and  s3(Ar). 

RIS-3.  e(Ar)  is  independent  of  s2(A;),  the  symbol  compo¬ 
nents  spanned  by  the  feedback  filter. 
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Although  not  defined  in  terms  of  e(k ),  RIS-1  was  used  in  [1] 
and  RIS-2  was  used  in  [24].  RIS-3  was  used  in  [4]  for  2-PSK, 
and  amounts  to  lumping  together  in  the  same  state  all  feedback 
error  sequences  e(A;)  0  s2(A;)  corresponding  to  the  same  e(k) 
and  assuming  they  have  equal  probability  of  occurring.  In  ef¬ 
fect,  because  the  distribution  of  the  error  sequences  is  unknow  n, 
we  make  the  approximation  that  they  are  uniformly  distributed 
This  is  the  key  assumption  in  this  paper  and  is  a  generalization 
of  [4]  to  M-PSK.  An  approximate  Markov  model  of  the  DFE 
follows  from  RIS-3  in  which  each  of  the  (2 M)Ad  states  is  de¬ 
fined  by  a  realization  of  e(k). 

To  compute  the  state  transition  probabilities  of  the  Markov 
matrix  P,  we  condition  y(k)  in  (19)  on  e(k)  representing  eaeh 
state.  The  transition  probabilities  from  the  pth  state  are  given  by 

Pqp  =  Pr(j/(k)  e  Sq\ep}  s(k)  =  1),  1  <  q  <  2 M  (20) 

where  ep  is  the  multiplicative  error  sequence  of  the  pth  state, 
and  Sq  is  the  soft  decision  region  in  the  eomplex  plane  corre¬ 
sponding  to  the  transmitted  symbol  s(k)  =  1  and  the  multi¬ 
plicative  feedback  error  eq.  We  construct  the  matrix  such  that 
the  transition  probabilities  from  the  pth  state  reside  in  the  pth 
row  of  P  Also,  the  first  row  of  P  corresponds  to  the  error- free 
state  in  which  caeh  element  of  ©i  is  ei  (Fig.  4). 

Once  we  have  constructed  P,  we  compute  the  steady-state 
probability  vector  v  as 


ptv  =  v 


(21) 


TABLE  I 

Channel  Model  With  Path  Powers  and  Normalized  Delays 


Path 

a'i  (dB) 

A,/T 

1 

-6.0 

-0.8 

2 

0.0 

0.0 

3 

-2.0 

1  8 

4 

-4.0 

2.7 

5 

-6.0 

3.6 

TABLE  11 

Static  Channel  Realizations 

Channel  X 

Channel  2 

r/T 

0 . 3837 

-0.4498 

Pi 

0.1916  +  iO . 3570 

0.0346  +  10.3766 

Pi 

0.4836  +  10.2146 

0.5033  ~  iO . 2058 

P3 

-0.3314  -  10 . 4589 

-0.8192  +  10.4337 

A 

-0.1165  -  iO .2191 

0.4052  +  iO.5288 

Pi 

0.5376  +  iO .3074 

-0.0024  +  10.1326 

DFE  is  given  by  Pse  =  1  —  p},  and  is  the  traditional  optimum 
measure  of  DFE  performance.  By  making  the  error-free  state 
absorbing  [8],  we  define  the  mean  burst  error  length  as  [25,  pp. 
43-52] 

Mbl  =  jtJni  -  (Nq  -  1)  (25) 

where  N  =  (I  —  G)-i,  and  ttq  is  the  initial  error  state  proba¬ 
bility  vector  that  we  define  as 


where  the  component  v (p)  represents  the  probability  of  being  in 
the  pth  state  [1].  Then  the  probability  of  M-PSK  symbol  error 
at  the  output  of  the  hard  decision  deviee  is  [16] 


(2A/pvn 

Te=  £  v{p)Pr(y(k)  i  s(k)  =  1)  (22) 

P—1 


*0  = 


1  -pE 


(26) 


Wc  include  the  term  N b  —  1  in  the  definition  of  /*bl  so  that 
the  minimum  mean  burst  error  length  is  1 .  Also,  the  variance  of 
the  burst  error  length  is  readily  computed  with  this  theory  and 
is  given  by  [8],  [25,  p  52] 


and  the  probability  of  bit  error  is 


4L=4(2N-I)N1-(4N1)2.  (27) 


(2M)nb  M 

Pa  =  log2  M  ^  ^ 

P=1  m=2 

xv(p)WmPi(y(k)  e  Hm \epy  s(k)  =  1)  (23) 


where  Hm  is  the  hard  decision  region  in  the  complex  plane  cor¬ 
responding  to  the  transmitted  symbol  s(k)  =  m,  and  Wm  is  the 
Hamming  weight  of  the  bits  assigned  to  symbol  m  with  respeet 
to  symbol  m  —  1.  Wc  calculate  the  probabilities  in  (20),  (22), 
and  (23)  by  averaging  over  the  noise  and  all  permutations  of  the 
symbol  vectors  sx (&),  s 2(fc),  and  s 3(k)  in  (19)  using  a  compu¬ 
tational  technique.  Details  are  in  the  Appendix. 

Next  we  partition  the  Markov  matrix  as 


P  = 


pi 


G. 


(24) 


If  we  set  7  =  1  and  0  =  0  (SL  =  H i),  then  the  probability 
of  symbol  error  assuming  error-free  feedback  for  the  standard 


Components  of  tto  represent  transition  probabilities  out  of 
the  error-free  state  weighted  by  the  probability  of  an  error  This 
definition  of  xq  is  convenient  for  high-order  modulation  and 
soft  decision  feedback  because  it  takes  into  account  that  certain 
errors  are  much  more  likely  than  others,  c.g.,  e2  and  eg  are  more 
likely  than  and  eg  for  4-PSK  in  Fig.  4. 

We  also  can  further  reduce  the  computational  complexity  of 
this  approximate  Markov  DFE  model  by  combining  states  con¬ 
taining  unlikely  errors  using 

RIS-4.  All  state  probabilities  in  (20)  that  transition  to  er¬ 
rors  eq  with  magnitude  that  exceed  some  threshold, 
are  lumped  into  the  worst-case  error  represented  by 
CAf  +  l- 

Consequently,  Markov  states  defined  by  that  differ  only  in 
elements  corresponding  to  these  unlikely  errors  are  incorpo¬ 
rated  into  an  aggregate  state  with  ca/+i  in  these  error  loca¬ 
tions.  For  instance,  for  4-PSK  we  might  use  only  the  errors 
{d,  e2,  63,  65,  C7,  eg}  in  Fig.  4  with  e.4  and  e6  incorporated 
into  65,  resulting  in  6Nb  Markov  states  rather  than  8Ab.  Wc 
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Fig.  5.  /*bl  versus  PE  for  Channel  1  with  4-PSK,  and  SNR  =  15  dB.  (a)  Over  soft  decision  angle  0  with  7  =  1.0.  (b)  Over  norm  constraint  7  with  6  =  10°. 


set  the  threshold  sueh  that  the  decrease  in  computational  com¬ 
plexity  is  worth  the  potential  loss  in  aeeuraey. 

V.  Numerical  Results 

To  demonstrate  the  utility  of  the  DFE  design  and  the  per¬ 
formance  analysis  approach,  we  present  results  for  two  statie 


channel  realizations  derived  from  the  ehannel  model  given  in 
Table  I.  The  fading  parameters  are  Rayleigh  distributed  with 
given  powers.  The  corresponding  normalized  path  delays  are 
also  shown.  The  equalizer  is  eoarse-synehronized  to  the  seeond 
path.  The  two  statie  channels  extracted  from  this  model  are 
given  in  Table  II,  and  present  different  degrees  of  challenge  to 
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(a) 


Fig.  6.  Theoretical  performance  curves  over  soft  decision  angle  (0)  and  norm  constraint  (7)  for  Channel  1  with  4-PSK.  (a)  SNR  =11  dB.  (b)  SNR  =  19  dB 


the  DFE.  The  normalized  timing  error  is  also  given.  The  DFE 
design  consists  of  a  fractionally  spaced  feedforward  filter  wF 
with  N  =  10  pre-  and  postcursor  taps,  and  a  feedback  filter  wn 
with  Nb  =  4  taps.  The  transmitted  square-root  raised  cosine 
pulse  has  30%  excess  bandwidth. 


A.  Channel  l :  4-PSK 

Wc  begin  by  equalizing  Channel  1  with  4-PSK  modulation 
and  cr|NR  =  15  dB.  Fig.  5(a)  displays  mean  burst  error  length 
ft bl  versus  probability  of  symbol  error  PE  as  the  soft  decision 
angle  0  varies  from  0  =  0°  to  0  =  30°  in  two-degree  incre- 
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Fig.  7.  Hbl  versus  PE  for  Channel  2  with  4-PSK,  and  SNR  =  12  dB.  (a)  Over  soft  decision  angle  0  with  7  =  1.0.  (b)  Over  norm  constraint  7  with  0  =  10°. 


ments.  No  norm  constraint  is  used  (7  =  1).  Also  shown  is  the 
probability  of  error  of  the  standard  unconstrained  DFE  given  no 
feedback  error  Fne-  Theoretical  and  simulation  results  taken 
over  2  x  108  transmitted  symbols  arc  displayed.  In  the  simula¬ 
tions,  error  bursts  are  determined  at  the  output  of  the  hard  de¬ 


cision  device  in  Fig.  2  and  are  defined  to  contain  no  more  than 
Nq  —  1  consecutive  correct  decisions.  However,  the  theoretical 
definition  of  /ibl  in  (25)  is  defined  with  respect  to  the  soft  deci¬ 
sion  device.  And  because  an  adjacent  soft  decision  error  (e2  and 
e8  in  Fig.  4)  may  not  necessarily  trigger  a  hard  decision  error. 
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Fig.  8.  Pk  versus  SNR  for  Channel  2  with  4-PSK. 


theoretical  and  measured  ^bl  tend  to  deviate  for  large  values 
of  0. 

Fig.  5(a)  reveals  that  the  traditional  DFE  using  no  soft  deci¬ 
sions  (0  =  0°),  suffers  degradation  in  Pe  with  respect  to  the 
error-free  feedback  case  P>je«  However,  as  we  increase  0,  an 
“L-curve”  [21  ]  develops  in  which  there  is  a  decrease  in  Pe  that 
is  accompanied  by  a  significant  decrease  in  /4bl  until  0  «  10° 
when  Pe  begins  to  increase.  Also  shown  is  the  considerable  im¬ 
provement  in  standard  deviation  of  the  burst  error  length  ctbl 
determined  theoretically  from  (27)  for  =  0°  and  0  =  10°. 

Next,  by  setting  0  =  10°,  we  begin  at  the  knee  of  the  L-curve 
in  Fig.  5(a)  and  plot  /*bl  versus  Pe  over  the  norm  constraint  7 
in  Fig.  5(b).  Wc  see  that  further  improvement  in  Pe  and  /*bl 
is  obtained  by  incorporating  the  norm  constraint  with  the  soft 
decision  device.  Again  an  L-curve  develops  in  which  both  Pe 
and  /ibl  decrease  until  7  »  0.70.  After  this  point,  any  decrease 
in  f  t bl  comes  at  the  expense  of  increasing  Pe.  At  the  operating 
point  7  =  0.70,  wc  have  decreased  Pe  approximately  by  half 
with  respect  to  the  standard  DFE,  and  have  come  closer  to  the 
optimum  Pne-  We  have  also  dramatically  decreased  ^bl  and 
<7bl-  For  the  purpose  of  comparison,  estimated  ^bl  and  <tbl 
derived  from  the  simulations  are  /7bl  =  7.1  and  <7bl  =  0.9 
for  the  parameter  set  (0  =  0°,  7  =  1.0),  and  fiui  =  2.4  and 
aBL  =  3.2  for  (0  =  10°,  7  =  0.70).  Wc  refer  to  the  knee  of 
this  L-curve  as  the  “optimum”  operating  point.  Also,  we  achieve 
the  same  result  if  we  generate  the  L-curves  first  over  7  and  then 
over  0. 

The  “optimum”  set  of  parameters  ( 0 ,  7)  generally  is  depen¬ 
dent  on  SNR.  To  demonstrate,  Fig  6  contains  three  dimen¬ 
sional  plots  of  theoretical  DFE  performance  for  Channel  1  with 
cr|NR  -  n  dB  and  cr|NR  =  19  dB  The  standard  deviation  l 


is  plotted  on  the  z-axis.  In  each  plot,  the  performance  curve  over 
0  with  7  =  1  is  combined  with  the  curve  over  7  with  0  =  10°. 
These  plots  reveal  that  although  wc  can  achieve  better  perfor¬ 
mance  at  the  low  SNR  using  a  different  parameter  set,  the  exact 
choice  of  0  and  7  is  not  critical  because  a  range  of  parameters 
produces  approximately  the  same  performance.  However,  the 
set  (0  =  10°,  7  =  0.70)  represents  a  good  compromise  be¬ 
tween  the  low  and  high  SNRs.  These  plots  also  demonstrate  the 
substantial  improvement  in  <tbl  that  can  be  obtained  using  the 
constrained  DFE  with  soft  decision  feedback,  particularly  at  the 
high  SNRs. 

B.  Channel  2:  4-PSK 

Fig.  7  contains  results  for  the  more  benign  Channel  2  with 
<r|NR  =  12  dB  Although  the  performance  gains  are  more 
modest,  wc  still  see  the  characteristic  L-curve  in  both  plots. 
Also,  it  is  interesting  that  (0  =  10°,  7  =  0.70)  is  close  to  the 
optimum  operating  region  of  the  DFE  for  both  channels  as  in¬ 
dicated  by  Figs.  5(b)  and  7(b).  This  was  observed  over  a  variety 
of  channel  realizations  for  4-PSK 

DFE  performance  is  given  as  a  function  of  SNR  for  the  two 
operating  regions  (0  =  0°,  7  =  1.0)  and  (0  =  10°,  7  =  0.70) 
in  Figs.  8  and  9  for  Channel  2  (similar  effects  arc  seen  for 
Channel  I).  Increasing  SNR  from  0  dB  has  the  effect  of  im¬ 
proving  Pe  and  /*bl  for  both  DFE  designs  because  in  this  re¬ 
gion  noise  is  the  major  cause  of  error  bursts  rather  than  decision 
error  feedback.  However,  after  approximately  1 1  dB,  the  mean 
burst  error  length  of  the  standard  DFE  in  Fig.  9  begins  to  in¬ 
crease  even  as  Pe  in  Fig.  8  decreases.  At  this  point  the  norm 
of  the  feedback  filter  increases  to  the  extent  that  decision  errors 
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Fig.  9.  ft bl  versus  SNR  For  Channel  2  with  4-PSK. 


TABLE  III 

Burst  Length  (BL)  Simulation  Statistics,  Channel  2,  4-PSK,  SNR  =  15  dB 


DFE  Type 

BL  -  1 

Number 

BL  «  2  BL 

of  Bursts 

>5  BL  >  10 

BL  >  15 

Maximum 

BL 

7=1,  0  =  0 

1356 

413 

517 

183 

61 

36 

II 

Qb 

II 

»—• 

O 

2945 

53 

29 

3 

0 

12 

begin  to  dominate  As  SNR  increases  beyond  1 1  dB,  initial  de¬ 
cision  errors  tend  to  cause  longer  and  longer  strings  of  errors. 
However,  bursts  occur  less  frequently  as  the  SNR  increases  to 
the  extent  that  the  longer  bursts  do  not  severely  degrade  the  av¬ 
erage  probability  of  error.  Conversely,  /*bl  of  the  constrained 
DFE  with  the  soft  decision  device  continues  to  decrease  as  SNR 
increases.  To  further  illustrate,  Table  III  contains  burst  length 
simulation  statistics  for  SNR  =  15  dB,  measured  over  10  trans¬ 
mitted  channel  symbols.  Although  the  constrained  DFE  exhibits 
about  twice  as  many  single  errors  as  the  standard  equalizer,  there 
is  a  significant  decrease  in  the  number  of  bursts  of  length  2  and 
greater.  Figs.  8  and  9  demonstrate  the  importance  of  ineluding 
burst  error  statistics  in  the  analysis  of  a  DFE. 

C.  Channel  2:  8-PSK 

This  example  demonstrates  DFE  performance  with  8-PSK 
modulation  for  Channel  2.  Fig.  10  displays  ^bl  versus  P e 
in  which  we  vary  0  and  7  with  a|NR  =  ^  dB.  To  min¬ 
imize  computations,  we  use  the  aggregate-state  approach 
described  by  RIS-4  at  the  end  of  Section  IV.  Adapting  the 
sequencing  from  4-PSK  in  Fig.  4  to  8-PSK,  we  use  the  errors 
{ci,  C2,  63,  e4,  69,  ei4,  ei5,  eie}  to  define  the  Markov  states, 
with  the  omitted  errors  lumped  into  eg.  We  then  have  8^® 


states,  which  is  the  same  number  used  to  generate  the  perfor¬ 
mance  results  for  4-PSK.  Comparing  Figs.  7  and  10  reveals 
that  we  need  to  increase  the  SNR  for  8-PSK  by  more  than  6 
dB  to  achieve  approximately  the  Pe  performance  of  4-PSK. 
Also,  the  burst  error  statistics  /^bl  and  <tbl  have  degraded 
significantly  for  the  standard  equalizer  with  (0  =  0°,  7  =  1) 
in  Fig.  10(a).  And  yet,  by  using  soft  decisions  and  invoking 
the  norm  constraint,  wc  sec  from  Fig.  10(b)  that  we  can  get 
/jbl  and  <tbl  values  comparable  to  those  given  in  Fig.  7(b) 
for  4-PSK  Also,  given  the  results  for  4-PSK,  it  is  intuitively 
satisfying  that  the  optimum  8-PSK  decision  angle  is  0  =  5°. 

D.  Channel  2:  4-PSK  With  Convolutional  Coding 

Although  the  increase  in  Pe  of  the  standard  DFE  with 
respect  to  Aje  may  appear  minor  in  Fig.  8,  the  associated  burst 
errors  have  significant  implications  for  error  correction  coding 
and  interleaving  [6].  This  is  particularly  true  when  operating 
at  high  SNR.  To  demonstrate,  we  incorporate  into  the  design 
a  simple  rate  1/2  convolutional  code  with  constraint  length  3 
and  free  distance  5  [26,  p.  466]  with  hard-decision  decoding. 
The  channel  symbols  arc  4-PSK  with  Gray  bit  mapping.  An 
( Ni ,  N{)  convolutional  interleaver  is  included  between  the 
output  of  the  encoder  and  the  4-PSK  channel  symbol  mapper, 
incurring  an  end-to-end  delay  of  Ni(Nj  —  1)  [27,  pp.  347-349]. 
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(b) 


Fig.  10.  /*bl  versus  PE  for  Channel  2  with  8-PSK,  and  SNR  =  18  dB  (a)  Over  soft  decision  angle  0  with  7  =  1.0.  (b)  Over  norm  constraint  7  with  0  =  5°. 


Fig  1 1  displays  the  probability  of  information  bit  error  Pi  at 
the  output  of  the  convolutional  decoder  as  a  function  of  SNR 
for  Channel  2  The  theoretical  upper  bound  on  the  error  of  the 
constrained  DFE  (0  =  10°,  7  =  0.70)  assuming  ideal  bit  inter¬ 
leaving  is  plotted.  This  bound  is  given  by  Pi  <  Aift(d). 


where  fid  =  (d  -  4)2d“5,  and  the  path  error  probability  /2(d) 
is  determined  using  the  channel  bit  error  probability  of  (23)  [26, 
pp.  463-466].  Also  included  are  simulation  results  for  the  stan¬ 
dard  DFE  and  the  constrained  DFE  with  convolutional  inter¬ 
leaving  and  Ni  =  11.  Although  the  constrained  DFE  mitigates 
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Fig.  ]  ] .  Pb  versus  SNR  for  Channel  2  with  rate  1/2  convolutional  code,  convolutional  interleaving,  and  4-PSK  channel  symbols. 


burst  errors,  it  does  not  eliminate  them.  And  because  convolu¬ 
tional  codes  are  sensitive  to  burst  errors,  an  interleaver  is  still 
required  to  achieve  the  modest  coding  gain  of  this  code.  How¬ 
ever,  for  1\  =  10-6,  approximately  2-dB  improvement  in  SNR 
is  obtained  with  the  constrained  DFE  versus  the  standard  DFE. 
To  achieve  this  same  performance,  iV*  must  double  for  the  stan¬ 
dard  DFE,  resulting  approximately  in  a  four-fold  increase  in 
end-to-end  delay. 


VI.  Conclusion 

We  have  demonstrated  that  a  DFE  design  that  combines  dis¬ 
crete  soft  decision  feedback  with  a  norm  constraint  on  the  feed¬ 
back  filter  is  an  effective  method  of  controlling  error  propaga¬ 
tion  with  minor  computational  cost.  We  analyzed  this  structure 
using  an  approximate  Markov  model  that  predicts  DFE  perfor¬ 
mance  in  terms  of  probability  of  symbol  error  as  well  as  burst 
decision  error  statistics  and  applied  it  to  M-PSK  modulation 

Other  soft  decision  implementations  than  the  simple  one 
proposed  here  for  M-PSK  may  offer  further  improvements 
in  performance.  For  instance,  we  can  quantize  the  decision 
space  more  finely  or  include  an  erasure  region  as  in  [16].  The 
applicability  of  the  analysis  approach  to  quadrature  amplitude 
modulation  (QAM)  can  be  explored.  Also,  we  may  achieve 
better  performance  by  using  more  complicated  DFE  designs 
that  incorporate  interleaving  and  error  correction  coding  into 
the  feedback  path  [28],  [29].  Finally,  we  have  considered  only 
individual  static  channel  realizations  drawn  from  a  fading 
propagation  environment.  Implementation  issues  related  to  an 
adaptive  DFE  operating  in  a  time-varying  environment  and 
the  effect  of  feedback  errors  on  tracking  performance  warrant 
further  analysis. 


Appendix 

Probability  Computations 


We  calculate  the  probabilities  in  (20),  (22),  and  (23)  using  a 
computationally  efficient  technique.  It  can  be  shown  that  [23] 

Pr(y(&)  €  D\ep .  s(k)  =  1) 

1  r  P  ~ 

~  ^  ^  Kd(vIi,  u>/q)/»isi(wZi,  wZq) 

°  iQ=-ri,=-r 

X  /lN(wZ i,  w/Q)eiMl^{ao}+u.-/QIm{a„})  +  (  (2g) 


where 

region  in  the  complex  plane  corresponding 
to  either  hard  or  soft  M-PSK  decisions; 
Fourier  transform  of  the  wedge-shaped  re¬ 
gion  D  truncated  at  T0j2  (see  Fig.  3  for 
D  =  50; 

characteristic  function  of  the  A/-PSK  ISI; 
characteristic  function  of  the  noise  process 
Wp  n(/c),  and  c o  =  2i r/T0. 

The  term  c  is  the  numerical  error  associated  with  this  method 
and  can  be  chosen  arbitrarily  small  with  the  proper  choice  of 
the  parameters  T0  and  P  [23]. 

The  characteristic  function  of  the  ISI  can  be  decomposed  ac¬ 
cording  to  the  symbol  vectors  si(A;),  S2(A;),  S3(/c),  and  the  as¬ 
sociated  coefficients  from  (19)  as 

d  Nb 

hisi(wlu  u)Iq)  =  ljIq)  /*2,  u>lQ) 

j= i 

d— A?q 

X  n  h3,iW,  u/Zq).  (29) 

i 


D 

u>Iq) 

ujIq) 

ui/q) 
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Each  component  of  the  1ST  corresponding  to  Si(k)  and  S3(A;)  is 
given  by 

hn,j(t*>lu  uIq) 

2  M/2  (  (2n  \ 

=  jg^co4w(Re{ani_,}/I+Im{a,l?_,}/y)cos(  — (m-  1)J 

m=l  '  '  / 

( 

+  ui{Re{anj}lQ  -  Im{aniJ}/i)  sin  I  _  1) 

(30) 

where  Re{a,l?J  }  and  Im{on  j}  are  the  real  and  imaginary  parts 
of  the  jth  clement  of  an  for  n  =  1,  3.  The  remaining  error 
feedback  terms  are  given  by 


/t2,j(uj/i,  w/q) 

2  M/2  ( 

=  ~m  ff  co?  (  w(Ile{o2,jer  j}h  +  Im{a2  jer,j}/y) 

m  —  1  ' 

Xcos(S(m_1)) 

+  w(Re{n2jer)^}/Q  -  hn{a2:jCp^}h) 

X  sin  (^(m  -  i)))  (31) 

where  epj  is  the  ^th  clement  of  the  multiplicative  error  vector 
ep  corresponding  to  thepth  state.  The  characteristic  function  of 
the  circular  Gaussian  noise  is  given  by 

/iN(wi,,  ulQ)  =  (32) 


where  2cr2  is  the  variance  of  the  complex  process  Wp  n (k). 

The  Fourier  transform  of  a  truncated  wedge  that  is  symmetric 
about  the  real  axis  is  given  by 


u>Iq) 


(— t/1  ^sine  Q  (ii  +  Iq  tan(</.))j 

X  e~inlQ  tan(0)/2 

-  sine  ($(/,-  Iq  tan (</>)))  j*1*  , 


/  Q^0 


K*{ 0, 


7T^ 

0)  =  tan(</>)^r 
u r 

^  _  i4tan(0)(-i)'' 
’  “  w2/2 


h^0 


where  0  is  the  angle  of  the  side  of  the  wedge  with  respect  to  the 
real  axis  and  sinc(x)  =  sin(7r:r)/7nr. 

For  example,  to  determine  the  components  of  (22)  for  4-PSK, 
we  calculate 


Pr(y(&)  $  Hi\ep,  s(k)  =  1) 

=  l-Pi(y(k)eH\\ep,s{k)  =  l)  (34) 

using  (28)  and  (33)  with  0  =  7t/4.  The  Markov  transition  prob¬ 
abilities  of  (20)  arc  determined  similarly  using  (33)  and  the  re¬ 
lation 

Ksq(uli,  u/Zq)  =  cos(<f>q)  +  ZQsin(0g)), 

u;(Zq  oos(0^)  -  l\  si n(0, )))  (35) 


where  <j>q  is  the  angle  of  rotation  of  the  soft  decision  region 
Sq  with  respect  to  the  real  axis.  For  the  4-PSK  soft  decisions 
defined  in  Fig.  3,  we  use  the  wedge  angles  <j>  €  {0/2,  7r/4  — 
0/2}  in  (33). 

References 

[1]  P  Monscn,  ‘Adaptive  equalization  of  ihc  slow  fading  channel,”  IEEE 
Trans.  Commun .,  vol.  COM-22,  pp.  1064-1075,  Aug.  1974. 

[2]  S.  U.  H  Qurcshi,  “Adaptive  equalization,”  Proc.  IEEE ,  vol.  73,  pp. 
1349-1387,  Sept.  1985 

[3]  P.  Monscn,  “MMSE  equalization  of  inlcrfcrcncc  on  fading  diversity 
channels,”  IEEE  Trans.  Commun .,  vol.  COM-32,  pp.  5-12,  Jan.  1984. 

[4]  L.  Li  and  L.  B.  Milstcin,  “Rejection  of  CW  interference  in  QPSK 
systems  using  decision-feedback  filters,  ’  IEEE  Trans.  Commun .,  vol. 
Com-31,  pp.  473-483,  Apr.  1983. 

[5]  K.  A.  Hamicd  and  G.  L.  Sliibcr,  “Performance  of  trcllis-codcd  modu¬ 
lation  for  equalized  muliipath  fading  1S1  channels,”  IEEE  Trans  Veh. 
Technol .,  vol.  44,  pp.  50-58,  Feb.  1995. 

[6]  D  J.  Cosicllo,  P  R.  Chcvillat,  and  L.  C.  Perez,  “Channel  coding  consid¬ 
erations  for  wireless  LAN’s,”  in  Proc.  IEEE  Int.  Symp.  Inform  Theory , 
1997,  p.  462. 

[7]  P.  Monscn,  “Theoretical  and  measured  performance  of  a  DFE  modem 
on  a  fading  multipath  channel,”  IEEE  Trans.  Commun.,  vol.  COM-25, 
pp.  1144-1153,  Oct.  1977. 

[8]  R  A.  Kennedy  and  B.  D.  O.  Anderson,  “Recovery  times  of  decision 
feedback  equalizers  on  noiseless  channels,”  IEEE  Trans.  Commun.,  vol. 
COM-35,  pp.  1012  1021,  Oct.  1987. 

[9]  D  L,  Duttwcilcr,  J.  E.  Mazo,  and  D  G.  Mcsscrschmitl,  “An  upper  bound 
on  the  error  probability  in  decision-feedback  equalization,”  IEEE  Trans. 
Inform.  Theory ,  vol.  IT-20,  pp.  490^-497,  July  1974. 

[10]  N.  C.  Beaulieu,  “Bounds  on  recovery  times  of  decision  feedback  equal¬ 
izers,”  IEEE  Trans.  Commun.,  vol.  42,  pp.  2786-  2794,  Oct.  1994. 

[11]  S  A.  Altckarand  N.  C.  Beaulieu,  “Upper  bounds  to  the  error  probability 
of  decision  feedback  equal  ization,”/£££  Trans.  Inform.  Theory,  vol.  39, 
pp.  145-156,  Jan  1993. 

[12]  W.  W.  Choy  and  N.  C.  Beaulieu,  “Improved  bounds  for  error  recovery 
times  of  decision  feedback  equalization,”  IEEE  Trans.  Inform.  Theory, 
vol.  43,  pp.  890-902,  May  1997. 

[13]  P.  S.  Bcdnarz  and  J  M.  Cioffi,  “Decision  feedback  equalization  for 
channels  with  error  correcting  capabilities,”  in  Proc  IEEE  Int  Conf. 
Commun.  (ICC  ’97),  1997,  pp.  1607-1612. 

[14]  G.  Mathew,  B.  Farhang-Boroujcny,  R  W.  Wood,  and  L.  Bin,  “Con¬ 
strained  equalizer  design  for  MDFE  detection  on  the  magnetic  recording 
channel,”  in  Proc.  IEEE  Global  Telecommun.  Conf.  (GLOBECOM  97), 
1997,  pp.  1258  -1262. 

[15]  R  A.  Kennedy  and  Z.  Ding,  “Quantizer  design  and  optimization  in  de¬ 
cision  feedback  equalization,”  in  Int.  Symp.  Signal  Processing  and  its 
Applications  (ISSPA-96),  1996,  pp.  614-617. 

[16]  M.  Chiani,  “Introducing  erasures  in  decision- feed  back  equalization  to 
reduce  error  propagation,”  IEEE  Trans.  Commun.,  vol.  45,  no.  7,  pp 
757  760,  July  1997. 

[17]  Y.  II  Kim  and  S.  Shamsundcr,  “Adaptive  algorithms  for  channel  equal¬ 
ization  with  soft  decision  feedback,”  IEEE  J  Select.  Areas  Commun., 
vol.  16,  pp.  1660-1669,  Dec.  1998. 

[18]  S.  L.  Ariyavisitakul  and  Y.  Li,  “Joint  coding  and  decision  feedback 
equalization  for  broadband  wireless  channels,”  IEEE  J.  Select.  Area 9 
Commun.,  vol  16,  pp  1670-1678,  Dec.  1998. 

[19]  J.  R.  Trcichlcr,  1.  Fijalkow,  and  C.  R.  Johnson,  “Fractionally  spaced 
equalizers,”  IEEE  Signal  Processing  Mag.,  vol.  13,  pp.  65-81,  May 
1996. 

[20]  P.  Balaban  and  J.  Salz,  “Optimum  diversity  combining  and  equaliza¬ 
tion  in  digital  data  transmission  with  applications  to  cellular  mobile 
radio — Part  1:  Theoretical  considerations,”  IEEE  Trans.  Commun.,  vol. 
40,  pp.  885  894,  May  1992. 

[21]  P.  C.  Hansen,  “Analysis  of  discrete  ill-posed  problems  by  means  of  the 
L-curvc,”  SIAM  Rev.,  vol.  34,  no.  4,  pp.  561  580,  Dec.  1992. 

[22]  S.  Haykin,  Adaptive  Filter  Theory,  3rd  cd.  Englewood  Cliffs,  NJ:  Pren¬ 
tice-Hall,  1996. 

[23]  M.  Reuter,  “Numerically  efficient  Fourier-based  technique  for  calcu¬ 
lating  error  probabilities  with  intersymbol  interference,”  IEEE  Trans 
Commun.,  vol.  45,  pp.  629-632,  June  1997. 

[24]  J  E.  Smec  and  N.  C.  Beaulieu,  “Error-rate  evaluation  of  linear  equaliza¬ 
tion  and  decision  feedback  equalization  with  error  propagation,”  IEEE 
Trans.  Commun.,  vol.  46,  pp.  656-665,  May  1998. 

[25]  J.  G.  Kcmcny  and  J.  L.  Snell,  Finite  Markov  Chains.  Princeton,  NJ: 
Van  Nostrand,  1960. 


REUTER  etal. :  MITIGATING  ERROR  PROPAGATION  EFFECTS  IN  A  DFE 


2041 


[26]  J.  Proakis,  Digital  Communications ,  2nd  cd:  McGraw-Hill,  1989. 

[27]  G.  C.  Clark  Jr.  and  J.  B.  Cain,  Error-Correction  Coding  for  Digital  Com¬ 
munications.  New  York:  Plenum,  1981. 

[28]  M.  V.  Eyuboglu,  “Detection  of  coded  modulation  signals  on  linear,  se¬ 
verely  distorted  channels  using  decision-feedback  noise  prediction  with 
interleaving,”  IEEE  Trans.  Commun.y  vol.  36,  pp.  401-409,  Apr  1988. 

[29]  K.  Zhou,  J.  G.  Proakis,  and  F.  Ling,  “Decision-feedback  equalization 
of  time -dispersive  channels  with  coded  modulation,”  IEEE  Trans. 
Commun.y  vol.  38,  pp.  18-24,  Jan.  1990. 


Michael  Reuter  (S’82  M’86)  received  the  B.S  and  M.S.  degrees  in  electrical 
engineering  from  the  University  of  Illinois,  Urbana-Champaign,  in  1984  and 
1986,  respectively,  and  the  Ph  D.  degree  in  electrical  and  computer  engineering 
from  the  University  of  California,  San  Diego,  in  2000. 

Since  1987  he  has  been  with  the  Space  and  Naval  Warfare  Systems  Center, 
San  Diego.  His  current  research  interests  arc  in  adaptive  and  statistical  signal 
processing  applied  to  problems  in  wireless  communications. 


Jeffery  C.  Allen  (M’94)  received  the  Ph.D.  degree  in 
mathematics  from  the  University  of  California,  San 
Diego,  in  1988. 

Since  1989  he  has  been  a  scientist  with  the  Space 
and  Naval  Warfare  Systems  Center,  San  Diego.  His 
current  work  is  focused  on  H***  applications  to  elec¬ 
trical  circuits,  stochastic  sea  surfaces  for  RF  propa¬ 
gation,  and  wireless  channel  modeling. 


James  R.  Zcidlcr  (M’76-SM’84-F’94)  received  the 
Ph.D.  degree  in  physics  from  the  University  of  Ne¬ 
braska,  Lincoln,  in  1972. 

Since  1974  he  has  been  a  scientist  with  the  Space 
and  Naval  Warfare  Systems  Center,  San  Diego,  CA. 
He  has  also  been  an  Adjunct  Professor  of  electrical 
and  computer  engineering  with  t  the  University 
of  California,  San  Diego,  since  1988.  During  this 
period,  he  conducted  research  on  communica¬ 
tions  systems,  sonar  and  communications  signal 
processing,  communications  signals  exploitation, 
undersea  surveillance,  array  processing,  underwater  acoustic  communications, 
infrared  image  processing,  and  electronic  devices.  He  was  also  a  technical 
advisor  with  the  Office  of  the  Assistant  Secretary  of  the  Navy  (Research, 
Engineering  and  Systems),  Washington,  DC  from  1983  to  1984. 

Dr.  Zcidlcr  was  an  Associate  Editor  of  the  IEEE  Transactions  on  Signal 
Processing  from  1991  to  1994  and  is  a  member  of  the  Editorial  Board  of  the 
Journal  of  the  Franklin  Institute.  He  received  an  award  for  the  Best  Unclassified 
Paper  at  the  IEEE  Military  Communications  Conference  in  1995  and  the  Navy 
Meritorious  Civilian  Service  Award  in  1991. 


Richard  C.  North  (S’84-M’90)  received  the  B.S.  degree  from  the  University  of 
California,  Davis,  in  1984,  the  M  Eng.  degree  from  Cornell  University,  Ithaca, 
NY,  in  1985,  and  the  Ph.D.  degree  from  the  University  of  California,  San  Diego, 
in  1991,  all  in  electrical  engineering. 

Since  1992,  he  has  been  with  the  Space  and  Naval  Warfare  Systems  Center, 
San  Diego.  From  1992  to  1994,  he  worked  on  the  research  and  development  of 
combined  space-time  adaptive  filtering  algorithms  for  advanced  sonar  and  com¬ 
munication  applications.  From  1994  to  1999,  he  was  the  chief  engineer  on  the 
High  Data  Rate  Linc-of-Sight  maritime  communications  project.  During  this 
time,  he  was  in  charge  of  the  development  and  demonstration  of  LOS  commu¬ 
nications  equipment  between  ships  and  ship-to-shorc  at  data  rates  up  to  sev¬ 
eral  Mbps.  Since  1998,  he  has  been  the  technical  lead  on  the  Digital  Modular 
Radio  (DMR)  program.  DMR  is  a  software  programmable,  multiband,  multi¬ 
mode  radio  designed  for  surface  and  subsurface  Naval  platforms. 


